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We point out a possible generation mechanism of non-Gaussian bubbles in the sky due to bubble 
nucleation in the early universe. We consider a curvaton scenario for inflation and assume that 
the curvaton field <j), whose energy density is subdominant during inflation but which is responsible 
for the curvature perturbation of the universe, is coupled to another field a which undergoes false 
vacuum decay through quantum tunneling. For this model, we compute the skewness of the curvaton 
fluctuations due to its interaction with a during tunneling, that is, on the background of an instanton 
solution that describes false vacuum decay. We find that the resulting skewness of the curvaton 
can become large in the spacetime region inside the bubble. We then compute the corresponding 
skewness in the statistical distribution of the cosmic microwave background (CMB) temperature 
fluctuations. We find a non- vanishing skewness in a bubble-shaped region in the sky. It can be large 
enough to be detected in the near future, and if detected it will bring us invaluable information 
about the physics in the early universe. 

PACS numbers: 98.80.Cq, 04.62,+v, 98.80.Es 



I. INTRODUCTION 

Inflation, a stage of accelerated expansion in the very 
early universe, is now widely accepted as part of the stan- 
dard evolutionary scenario of the universe. On the other 
hand, many models of inflation have been proposed but 
we are still far from being able to narrow down the pos- 
sible models sufficiently. Among those many models of 
inflation, much attention has been paid recently to the 
ones based on string theory which is considered to 
be a promising candidate for the ultimate unified theory. 
In particular, it is of great interest if the string theory 
landscape 0], in which there are many local minima, or 
false vacua, and the universe jumps from one minimum 
to another by quantum tunneling, can be observationally 
tested 0, [if . Quantum tunneling of a scalar field with 
gravity is usually treated with the Coleman-De Luccia 
(CDL) instanton method [f|, where the evolution of a 
scalar field is described with an 0(4)-symmetric instan- 
ton, which is a solution of the Euclidean equations of mo- 
tion. Motivated by the string landscape, inflation models 
with multi-scalar fields and/or with tunneling are now 
keenly studied. Extension of the CDL instanton method 
to a multi-scalar field system is already discussed in @. 
In those studies, however, many inflation models have 
been proposed, and now it is important to distinguish 
those models by observation. 

Observations of the power spectrum of CMB temper- 
ature anisotropics have convinced us of the existence of 
an inflationary stage in the very early universe. On top 
of that, if a non-Gaussian feature such as skewness or 
bispectrum is detected in the CMB anisotropy, it will 
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have a strong impact on the physics of the early uni- 
verse Q. In particular, since any single- field slow- roll in- 
flation model with canonical kinetic term predicts almost 
Gaussian fluctuations 0, any non-zero non-Gaussianity 
will exclude all these models. Observations of non- 
Gaussianity use templates, such as local type [9j.[l9j 
equilateral type , and orthogonal type 11 1 , in order to 
increase the statistical significance. However, since all of 
these templates assume statistical isotropy, there may be 
anisotropic non-Gaussianities which may not have been 
detected by these templates. 

In this letter, we study a multi-field model in which 
a nonlinear interaction between two scalar fields, one 
of which being responsible for the curvature perturba- 
tion of the universe (that is, for the formation of the 
large scale structure) and the other for quantum tun- 
neling via a CDL instanton, induces an anisotropic non- 
Gaussianity. To be specific, we introduce an inflaton field 
<I> that realizes slow-roll inflation, a tunneling field a that 
governs the tunneling dynamics, and a curvaton field <fi 
that contributes to the curvature perturbation of the uni- 
verse mm. 

The inflaton dominates the energy density of the uni- 
verse during inflation but rapidly decays to radiation af- 
ter the end of inflation, On the other hand, the energy 
density of the curvaton field is negligible during inflation 
but its decay is delayed after inflation so that it gradually 
begins to dominate the universe. 

Approximating the universe during inflation by an ex- 
act de Sitter spacetime, the inflaton behaves as a cos- 
mic clock and determines an appropriate time-slicing, 
namely, a spatially flat time-slicing of the de Sitter space- 
time. In this setup, assuming that the energy scale as- 
sociated with the tunneling field is much smaller than 
the energy scale of inflation, the bubble nucleation can 
be well described by a single-field CDL instanton with 
no backreaction to the geometry, that is, on the exact de 
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FIG. 1: Penrose-like diagram of bubble nucleating universe, 
where a half of Euclidean region (bottom) and Lorentzian 
region after bubble nucleation (top) are shown The blue dot- 
dashed lines are surfaces of constant cosmic time, the green 
solid lines those of uniform tunneling field, and the red solid 
line the location of the bubble wall. 



Sitter background 0, [1] • 

The curvaton <j> is affected by the background bubble 
through a coupling with the tunneling field a. For sim- 
plicity and definiteness, we consider a potential of the 
form, Vint(c, 4>) = A(er)0 3 . Assuming that Vi n t(cr, <f>) is 
non-vanishing only at or inside the bubble wall, we ex- 
pect that 4> may have a spatially localized, bubble-shaped 
non-Gaussianity due to the background bubble-shaped 
configuration of a. This leads to an anisotropic, bubble- 
shaped skewness of the CMB temperature anisotropy. 

This paper is organized as follows. We first illustrate 
the background spacetime and the configuration of the 
bubble. Next, we briefly review a useful formalism for 
computing the equal-time iV-point functions, the tunnel- 
ing in-in formalism, and calculate the skewness of the 
curvaton fluctuations. Then we demonstrate that a sky 
map of an anisotropic non-Gaussian parameter /nl in our 
model. Finally, we end with a few concluding remarks. 



II. BACKGROUND EVOLUTION 

Let us start from a brief description of de Sitter space- 
time 0, EH • A Lorentzian 4-dimensional de Sitter space- 
time is 0(4,l)-symmetric, which is obtained by analyti- 
cal continuation of an 0(5)-symmetric 4-dimensional Eu- 
clidean sphere, as illustrated in Fig. [T] Among various 
coordinatization of the spacetime, the uniform tunnel- 
ing field slicing is most appropriate to see 0(4)(0(3,1))- 
symmctry of the CDL instanton, which corresponds to 
the time-slicing inside the bubble that describes a spa- 
tially homogeneous and isotropic open universe. For 
brevity, let us call it open slicing. Open slicing in the 
C-rcgion in Fig. Q] is a time-like slicing, an the metric in 



the C-region may be expressed as 
1 



ds z = 



H 2 



d X 2 



sin 2 x l 



dr z + cosh t dQi 



(1) 



and dVl\ is the 



where — oo < r < oo and < \ — n 
metric on the unit 2-sphere. Theses coordinates can be 
extended to the other regions of the spacetime by ana- 
lytical continuation: 



t = j'r + 27r/2 = ?'l — iir/2 
X = = -Hl = *e , 



(2) 



where (iR,rR.), and (tE,^) are the coordinates 

for the R-, L-, and E-regions, respectively, in Fig. [T] 
The metric for the flat time-slicing is 



ds 2 = -dt 2 



\t) (dr 2 



l dVLi 



a(t) = H~ 



(3) 



As mentioned in the introduction, the t = const, slices 
are on which the inflaton is uniform, and give the cosmic 
rest frame of our universe. 

Let us briefly describe the bubble configuration of a 
under the thin- wall approximation. The 0(4)-symmetric 
CDL instanton for a can be described as a function of Xi 
which we denote by a(x) hereafter. We denote the wall 
radius by i?w- Hence 
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for 



< X < HR W , 
HR W <x<k . 



(4) 



Though a is homogeneous in the R- or L-regions, the 
physical bubble radius £w{t) on flat slices increases as 
time goes on. We take the origin r = t = of 
the flat slicing metric <|3j> to be the center of the bub- 
ble at the time of nucleation. Then the relations be- 
tween the coordinates in the metrics ([l) and ([3]) are 



cosx = cosh-Hi — (l/2)e r 2 and sinxcoshr = e 
which gives £w(t) = a(t)\/l~ 
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As seen from this expression, a bubble once nucleated 
with radius £w(0) (~ Rw for i?i?w "C 1) expands 
to the Hubble horizon scale within one or two e-folds 
of time, At ~ and then expands comovingly as 

e w (t) ~ a(t). 

It should be noted that models with multiple nucle- 
ation is possible. If those bubbles do not interact each 
other, we can take into account the effect of all bubbles 
by summing up the effect of each bubble. Hereafter, we 
consider a model with a single bubble for simplicity. 



III. NON-GAUSSIANITY GENERATION 
A. Effective Action on Instanton Background 

Here we calculate the skewness in the curvaton fluctu- 
ations on the single bubble background. We consider the 
Lagrangian for the curvaton <j> as 

/ 1 TO 2 \ 

-9 ( ^ V d^d^ + — 2 + M„t(a, 4>) , (5) 
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where m 2 <C H 2 . We expand the above by setting <fi = 
4>q + 54>, where <j>o is a homogeneous classical part and 5(f) 
is a quantum fluctuation. We assume that 4>o is approxi- 
mately constant during inflation, and we concentrate on 
the evolution of 5(f). Further, for simplicity, we assume 
Vi n t (c, </>) is non- vanishing only on the wall. Hence we 
approximate it as VS n t(<7, 4>) = \<f> 3 H S(x — HRyj) with 
XH = JT dxX(a(x)) ■ Then the Lagrangian for 5<f> is 
given as Cg^ = £ + Cj , where Cq is free-part and Cj is 
interaction-part, 



( \aTW 



>d„ 



+ -(m 



2 ' 5m 2 ) 



d = -^jXH5(x- HR w )5<t> s 
where 5m 2 = 6XH(j) S(x — HR W ). 



(6) 



B. Quantum Field Theory on the Instanton 
Background 



To calculate the quantum fluctuations on the instanton 
background, we extend the in-in formalism to Euclidean 
spacetime, which we call the tunneling in-in formalism. 
This formalism is based on the WKB analysis of a tun- 
neling wave function for free theory [TEI, to the case 
with nonlinear interactions. A detailed derivation will be 
given elsewhere fl7j.(2Cl] 

The tunneling in-in formalism tells us that the iV-point 
function of 5<f> is given by 
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FIG. 2: Same figure as Fig. [T] but with all domains of inte- 
gration given in eq. J7J. 



(0\5(j)(x)5(j)(x')\0) using Wick's theorem. We note that 
G{x,x') is not the correlation function for the Bunch- 
Davies vacuum due to the non-trivial nature of £ . It 
may be obtained by studying the "evolution" of the mode 
functions in the Euclidean space 0] . 

While G(x,x / ) is a single- valued function when x and 
x' are in space-like separation, when they are in time- 
like separation, T or T in the expression singles out the 
Feynman or anti-Fcynman propagator, respectively. A 
branch of G(x, x'), when x' is in the Euclidean region, 
is determined by analyticity on x' along C. It may be 
noted that by this way of choosing branches the expres- 
sion in eq. ([7]) is equivalent to that obtained by analytical 
continuation of Euclidean quantum field theory [l~8j , and 
0(4)(0(3,1)) symmetry of the result is guaranteed. 
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(7) 

where the tunneling in-in path C and t = const, sur- 
faces £t are as shown in Fig. [2] The first half of C (the 
arrowed green line in Fig. [5]) goes from one end of the Eu- 
clidean region to future infinity in the Lorentzian region 
through the bubble nucleation surface. The second half 
of C (the arrowed blue line in Fig. [5]) goes back from fu- 
ture infinity through the nucleation surface to the other 
end of the Euclidean region. The slicing C x St covers 
the whole Euclidean region and the future half of the 
Lorentzian region twice. In the Lorentzian region any T, t 
is a Cauchy surface. The operator P in eq. ([7]) is the path- 
ordering operator. In the Lorentzian region, P reduces 
to the time-ordering and anti-time-ordering operators T 
and T on the first and second halves of C, respectively. It 
should be noted that eq. (UJ is independent of the choice 
of coordinates. 

Eq. ([7J is evaluated in the same way as in the 
usual perturbation theory. After expanding the inter- 
action part perturbatively, operators are transformed 
to products of the free correlation function G(x,x') = 



C. Skewness from Bubble Nucleation 

Now we evaluate the skewness (5(j) 3 (x)} by substituting 
in eq. © to the tunneling in-in formula in eq. (J7J. 
In the following calculations, we approximate G{x 1 x') by 
that of the Bunch-Davies vacuum, by neglecting correc- 
tions due to the non-trivial part of Co in eq. © • This ef- 
fect may affect the value of the skewness, but since it only 
induces a statistically homogeneous non-Gaussianity, we 
simply ignore it here. 

The free correlation function G(x,x') of the Bunch- 
Davies vacuum is given by 14 1 



G(x,x')-- 



167T COS 7TV 



1*1 



v 1 + Z(x,x') 



(8) 



where v = \/9/4 — m 2 /H 2 and Z(x,x') = 
cos(Hd(x,x')), where d(x,x') is the geodesic dis- 
tance between x and x' . We have Z > 1, Z = 1 and 
Z < 1, for a timelike (imaginary d), null (d = 0), and 
spacelike (real d) separation, respectively. As mentioned 
before, a branch of G(x,x') = G(Z(x,x')) should be 
specified when Z > 1 and it can be done by adding small 
imaginary part ±ie to Z as G(Z(x, x') ± ie) . 

We assume A < 1 to avoid possible strong coupling 
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FIG. 3: The structure of integrand and integration path in 
cq. (jQj . The wavy line is the branch cut, the dotted line is 
the original path C, while the solid line is a detoured path C' 
used for numerical evaluation. 



problems. To leading order in A, we obtain 

(Stfix) J> = -i\H~ 3 J dr' J d X ' J dfi' 

x sin 3 x' sinh 2 r' S ( X ' - HR W ) (G(x, x')f . (9) 



where 



C : 



7T 

-,--,0 



+ ie 



IT 

l 2 



(10) 



as shown in Fig. |3] Among the domains of integral, r' € 
(— i ir/2,+iir/2) and r' g (0 ± z e, oo) correspond to the 
E- and C-regions, respectively. Integration in the R- or 
L-regions is unnecessary because Cj vanishes there. A 
small imaginary part in r' € (0 ± ie, oo) gives a small 
imaginary part to Z(x, x') to select the Feynman or anti- 
Feynman propagator. 

Evaluation of eq. ([9]) is straightforward in principle, 
but it is significantly simplified if we use the 0(4)- 
symmetry of the Appoint functions in the tunneling in- 
in formalism. Just for illustration, let us consider the 
case when x is in the R-region. In this case, computing 
(S(j) 3 (x)) at = 0, where integration over x an d f2 is 
trivial, is enough since the value at any other point in 
the R-region can be known from the 0(4)-symmetry. 

After integrating over x and fi, the integration over r 
along C remains. Since the original path C passes near 
the singularity, it makes the integral difficult to evalu- 
ate numerically To avoid it, we deform the path C to 
another path C without crossing the branch cut or the 
poles of the integrand as in Fig. [21 The resulting (5cj) 3 (x)) 
is obtained in terms of the open-slice coordinates, but 
the expression in terms of the flat-slice coordinates can 
be easily found by the coordinate transformation. 

Fig. @] shows the skewness on flat slices at time < as a 
function of the radial coordinate r, (S(j> 3 (r,t)) / (5cj) 3 (0, t)), 
for Ht = 50 and ifi?w = 0.2-7T for several values of m/H. 
To see the dependence of the skewness on the model pa- 
rameters, we introduce a normalized skewness F(t) at the 
center of the bubble by 
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FIG. 4: (5(f> 3 (r,t)} normalized by the value at the center as a 
function of r at Ht = 50. 



F(t) 

o.oi 

io- 4 

10" 6 
10 

io-'°t lu- 
ff 





v Ht=50 


Ht=60 x 





(S(t) 3 {0,t)) = 4ttA sm 3 (HR w ) 



3H 2 
8n 2 m 2 



H 3 F(t). (11) 



FIG. 5: A normalized skewness as a function of time, F(t). 
The solid lines are the result of numerical evaluation and the 
dashed lines are analytical approximation, e~'a 7Ht . 



The result is shown in Fig. [5j From this, we see that 
F fa e 

We see that the skewness is large inside the bubble 
and it decreases as one moves away from the bubble, and 
its radial dependence is stronger for larger m/H. These 
features can be understood from the fact that G{x,x') 
is roughly given by (3iJ 4 /87r 2 m 2 )|l - Z\- m2 ^ 3H ^ for 
m/H <C 1 and Z is not very close to 1. The appar- 
ent divergent behavior at the bubble wall comes from 
the UV divergence in G(x,x'). But, this divergence dis- 
appears as soon as the finite resolution of observations or 
renormalization of the theory is taken into account. 



IV. NON-GAUSSIAN BUBBLES IN THE SKY 

Now we study signatures of our model in the CMB 
anisotropy. In our model, the evolution of the universe 
after the inflaton and the tunneling field decay, is exactly 
the same as the ordinary curvaton scenario hj. Here, we 
estimate the curvature perturbation after the curvaton 
decay by using the sudden-decay approximation fl3j . 

After inflation the curvaton starts to roll down the 
potential and undergoes damped oscillations when the 
Hubble parameter becomes smaller than m. The cur- 
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vaton energy density behaves like a pressurclcss matter 
during this stage of damped oscillations. At the time 
t = icurv when the curvatons decay, the energy density 
of the universe consists of that of radiation p r generated 
from the decay of the inflatons and that of the curvaton 
Prf, . The contribution of each component to the total cur- 
vature perturbation of the universe may be conveniently 
expressed in terms of £ r and C^, which are defined as the 
curvature perturbations on the uniform energy density 
slices of the radiation and the curvaton, respectively. 

It is known that each of them is separately conserved 
on superhorizon scales at t < t curv . As for C^, it may 
be evaluated by the energy density fluctuations on the 
uniform inflaton field slice at the end of inflation, i 
as 



3 I 2 d> 



(12) 



where t e is the time at the end of inflation. 

After the curvaton decay, the total curvature pertur- 
bation on the slices of uniform total energy density, £, is 
conserved on superhorizon scales. It is given by 



C = (1 - r^Cr + 



(13) 



where = 3/^/ (4p r + 3/ty)| tt!U1 . v . 

For the moment, we ignore the contribution from 
5(j) 2 /(j)Q and concentrate on the skewness generated by 
the nonlinear interaction with the tunneling field. Then 
we have 



8 3 W> 3 Mc)> 



27 



(14) 



where (5(f> 3 ) is given by Eq. ([9]). An important feature 
is that the skewness of the curvature perturbation de- 
pends on the position of the bubble due to the spatial 
dependence of (<5</> 3 ). 

To proceed, we focus only on the large angular scale 
CMB anisotropy for which the Sachs- Wolfe effect domi- 
nates. In this regime, the CMB temperature anisotropy 
can be written in terms of the curvature perturbation as 
(ST/T)(n) = (1/5) ((x + r*n,t*) , where x is the co- 
moving position of the observer measured from the bub- 
ble center, h = (sin 9 cos ip, sin 9 sin ip, cos 9) is the direc- 
tional cosine of the sky seen by the observer, is the 
comoving distance from the observer to the last scatter- 
ing surface and £* is the time at the last scattering sur- 
face. Then, defining the non-Gaussianity parameter /nl 
as (3/5)/ NL (n) = {C 3 {x + r*n, £*))/(C 2 ) 2 , we have 



hh(n) 



40r 3 (6cj) 3 (\xo+r,h\,U)) 
81(C 2 ) 2 $ 



(15) 



One explicitly sees that the dependence on the position 
of the bubble center breaks the statistical isotropy. 

For completeness, let us discuss the additional contri- 
butions to /nl which we have ignored. The contribution 
of Q r to non-Gaussianity is known to be very small as long 
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FIG. 6: Non-Gaussianity map in the CMB sky, for the model 
parameters given in the text. Colors correspond to the values 
of /nl. 



as the vacuum is in the Bunch-Davies vacuum |8( ■ A de- 
viation from the Bunch-Davies vacuum may give rise to 
a non-zero non-Gaussianity. It is expected to be statisti- 
cally homogeneous and isotropic but scale-dependent (in 
the Fourier space). It may be detected by the templates 
of the equilateral or orthogonal types. The contribution 
from the part of quadratic in 8<f>/ '4>o gives rise to a local 
type non-Gaussianity which is again statistically homo- 
geneous and isotropic, and which may be detected by the 
squeezed type templates. 

This contribution from the nonlinearity of S<f>/(j) in ^ 
may be estimated as follows. Roughly speaking, (<5</> 4 ) « 
(S(j) 2 ) 2 ~ (H/2t:) 4 for m/H <C 1, and the contribution to 
(C|) is about (4/9)(i//27T0 o ) 4 . Thus, the condition that 
/nl is dominated by the nonlinear interaction with the 
bubble is given as 



4 / H y 32tt . 3/rrn . / 3H 2 H 



This is satisfied in an example we compute below. 

For illustration, we plot /nl("-) in the CMB sky in 
Fig. |51 The parameters arc m/H = 0.3, i?i?w = 0.2-7T 
H/cf) = 0.001 A = 0.005, r = 0.1, |x | = r* = 2, 
Ht c = 50. For the variance (£ 2 ) we simply impose the 
observational result, (C 2 ) = A 2 C = 6.25 x 10~ 10 0, as- 

/.2\ ~ 



211 Since 



suming that it is dominated by Q r 
(H/2-k) 2 when m/H <C 1, this assumption is justified if 
(2r c i > /i) 2 (H/2'!T(j) ) 2 <C A 2 , which is marginally satisfied 
in the above example. 

The typical value of /nl at the center of the bubble 
when the parameters satisfy the above conditions is esti- 
mated as 



f (cen) 



3xlO~ 4 A7^ sin 3 (HR X 
Af exp((%) 2 Ht e ) 



m ) Go) ' ^ i6 ^ 



This gives /^ c ™' ~ 15 for the parameters given above and 
agrees with the result in Fig. [6] within the errors shown in 
Fig. [S] From this estimation, we find that the resultant 
/nl is rather sensitive to the values of m/H and H/4>q. 
It is enhanced for smaller m/H and larger H/(f>Q. 
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V. CONCLUSION 

In this paper, we calculated the skewness in the CMB 
temperature anisotropy in a model with bubble nucle- 
ation during inflation, motivated by the string theory 
landscape. We considered bubble nucleation in the cur- 
vaton scenario of inflation in which the curvaton vac- 
uum fluctuations are affected by the bubble nucleation 
through interaction with the tunneling field. The cal- 
culation was done by extending the in-in formalism to 
the instanton background [l7j]. We found that there can 
be spatially localized, bubble-shaped skewness which is 
large inside the bubble. 

As far as we know, bubble-shaped non-Gaussianitics 
have not been studied carefully yet in observation. So 
it seems interesting to look for such a non-Gaussianity 
already in the current observational data. As pointed 
out by Komatsu and Spergel Q, the bispectrum, cor- 
responding to the 3-point function, contains much more 
information than the single value of the skewness. Thus, 
analysis beyond skewness may improve the observability 
of bubble-shaped non-Gaussianities. We hope to come 
back to this issue in a future publication. 

Since the string theory landscape gives a strong moti- 
vation for inflation models with bubble nucleation, stud- 
ies of such inflation models may be regarded as testing 



string theory using the universe as a laboratory. In any 
case, if any signature of bubble nucleation during infla- 
tion is found in observation, it gives a huge impact on 
the physics of the early universe, including string theory. 

Finally, we note that in the evaluation of the skewness 
we neglected the effect of deviations from the Bunch- 
Davies vacuum. This may affect details of our result, 
though generic features are expected to remain the same. 
This effect can be evaluated by studying the evolution of 
mode function on the instanton background Q . We plan 
to come back to this issue in the near future. 
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